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Abstract 

The decay of a small homogeneous perturbation of the temperature of a dilute granular gas in 
the steady uniform shear flow state is investigated. Using kinetic theory based on the inelastic 
Boltzmann equation, a closed equation for the decay of the perturbation is derived. The equation 
involves the generalized shear viscosity of the gas in the time-dependent shear flow state, and 
therefore it predicts relevant rheological effects beyond the quasi-elastic limit. A good agreement 
is found when comparing the theory with molecular dynamics simulation results. Moreover, the 
Onsager postulate on the regression of fluctuations is fulfilled. 

PACS numbers: 45.70.-n,05.20.Dd, 05.60.-k,51.10.+y 



1 



I. INTRODUCTION 



The study of the linear response of an equilibrium system to a small perturbation either of 
the hydrodynamic fields or of an external field, is a formally exact method of non-equilibrium 
statistical mechanics, that has led to a significant advance in the theory of simple atomic 
liquids [l]. Nevertheless, little is known with the same degree of generality and exactness 
about the linear response of fluids that are in non-equilibrium steady states. Granular gases 
provide an appropriate context in which to address this as well as many other general ques- 
tions about irreversibility {2, s|. This is because granular gases are intrinsic non-equilibrium 
systems, but also because they exhibit steady states that are macroscopically simple, as 
compared with those of molecular fluids. 

The objective of this work is to provide a case in which the theory can be developed 
in a controlled way and to compare its predictions with molecular dynamics simulations. 
The linear response of a dilute granular gas in a reference homogeneous steady state with 
constant shear rate is studied, for arbitrary values of the latter. Although this state can not 
be implemented in real experiments, it has the advantage of being easily obtained by means 
of particle simulation methods. The considered perturbation is homogenous and affects the 
internal energy, or granular temperature, of the system but not the macroscopic velocity 
field. Rheological effects are known to be relevant in the reference state when the shear rate, 
_d in units of the inve.se of the average tnne between collisions, ,s nn.te ffl. These 
effects manifest themselves in the strong dependence of the components of the pressure 
tensor on the shear rate, beyond the proportionality relation assumed by the Navier-Stokes 
approximation. A question addressed here is how this rheology translates into the linear 
response of the system to the perturbation. 

A key ingredient in the linear response theory for molecular fluids is the Onsager postulate 
g asserting that the average regression of thermal fluctuations in equilibrium foflows the 
macroscopic linear law of relaxation. This property is now quite well understood but 
little is known about its validity in far from equilibrium states. This issue will be investigated 
here for a granular gas under steady shear flow. 

The remaining of this paper is organized as follows. In Sec. [Til the equations character- 
izing a dilute granular gas under uniform shear flow are presented, and the existence of a 
steady solution is indicated. A small homogeneous perturbation of the temperature of the 
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steady state is considered, and the linear equation describing its time evolution is derived. 
The scaling verified by the distribution function of the system is indicated and its conse- 
quences on several macroscopic properties of the relaxing system are discussed. This leads 
to a formal equation for the decay of the temperature. This equation is transformed into 
a closed hydrodynamic form in Sec. IHH by employing known properties of the reference 
state. Then, the theoretical prediction is compared with molecular dynamic results, and a 
good agreement is obtained, although some systematic discrepancy shows up for very strong 
inelasticity. Moreover, it is shown that the linear perturbation decays with the same rate as 
the temperature time- correlation function in the reference state, i.e. the Onsager postulate 
is verified. Finally, Sec. [TV] is used to present a short summary of the results obtained in 
the paper. 



II. DECAY OF A TEMPERATURE PERTURBATION 

The focus here is on the decay of a small homogeneous perturbation of the temperature in 
a granular gas in the steady uniform shear flow (USF) state. Although explicit calculations 
will be restricted to the low density limit, their extension to higher densities may be possible. 
Consider first a granular gas in the time-dependent USF state. At the macroscopic level, this 
state is characterized by a uniform number density n, a time-dependent uniform granular 
temperature T(t), and a constant velocity field u with linear profile 0, [j, 0, [9], 



^ j ®ij j 1 ^ij — Q'fiixfijyy (1) 



Ui\r r 
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where 6ij is the Kroneker delta function and a the constant shear rate. For this state, the 
mass conservation is identically verified, the conservation of momentum requires that the xy 
component of the pressure tensor P xy is uniform, and the energy balance equation reads [5] 

m = -* PM - m . (2) 

Here, d is the dimension, 2 or 3, of the system, and ((t) is the cooling rate due to the energy 
dissipation in collisions. The above relation admits a steady form, 

^P X y,s = -(sT s . (3) 

Equations ([2]) and ([3]) are exact consequences of the balance equations, but they become 
meaningful only after some constitutive relations for P xy and ((t) are given. These relations 



express the cooling rate and the pressure tensor as functionals of the hydrodynamic fields, 
n, u, and T. The particular case of a low density granular gas composed of smooth inelastic 
hard spheres (d = 3) or disks (d = 2) of mass m and diameter a will be considered. For 
states that are homogeneous in the Lagragian frame of reference for the velocity field given 
in Eq. (frj, the one particle distribution of the gas, f(V,t), where V(r) = v — u(r) is the 
peculiar vecity of the particle, obeys the inelastic Boltzmann equation 



^-*V y -^r)f(V,t) = J[V\f(t)], 



(4) 



where J 



V\f(t)] is the Boltzmann collison operator for smooth inelastic hard spheres or 
disks 10] . The shear stress P xy and the cooling rate ((t) in Eq. (j2J) are functionals of the 



one-particle distribution function, 



P xy (t) = m j dVV x V y f(V,t), 



at) 



(1 _ a 2 )7r (d-l)/2 mtr d-l 

4r (^) nT{t)d 



dV x / dV 2 V? 2 /(Vi,t)/(V 2 ,t), 



(5) 



(6) 



with V12 = V\ — "Vi- A sufficient condition for the existence of constitutive relations and 



hence of an hydrodynamic descri 
of the kind named normal 
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ption, is that the distribution function of the system be 



12j, and meaning that all the time dependence in it occurs 



through the temperature, i.e., it has the form 

f(V,t)=n[v (t)]- d X (c,a*), 



where 



is a thermal velocity, 



and 



V (t) EE 



a = 



2T(t) 
m 

V 

aX 



lV2 



(7) 

(8) 
(9) 
(10) 



with A 



[no 



being a quantity proportional to the mean free path of the gas. The 



distribution function of the steady USF state, is obtained by particularizing the above 
expression for the steady temperature, namely 



X(c,a*) -»■ x(c,a*), 



(11) 



'n., V "' 

The idea that for a wide class of initial conditions the steady USF state is reached, 
after a short initial transient period, through a hydrodynamic evolution of the system has 
been already used. In particular, the way in which the normal solution is approached has 



been investigated [13|, Il4j . Suppose that the system is in the steady USF state and then 
it is submitted at t = to an instantaneous and homogeneous small perturbation of the 
temperature (or internal energy), keeping constant the value of the shear rate a. It is assumed 
that the subsequent time evolution of the temperature is described by hydrodynamics, so 
that it will obey the linearization of Eq. ([2j) around the steady USF state. Such a linearization 
reads 

= ~Vd 6Pxy(t) ~ TM{t) ~ Cs5T(t) ' (13) 

where ST(t) = T(t) — T s , and SP xy and 5( are the linear deviations of the shear stress and 
the cooling rate, respectively, from their steady values. In the hydrodynamic regime, Eq. 
(JHJ) can be rewritten as 

P xy {t) = ±nmv 2 (t)P: y (a*), (14) 
with the reduced shear stress defined as 

P: y (a*) = 2 J dcc x c yX (c,a*). (15) 



Then, a simple calculation gives 



M*) 1 2 *f9P* y (a*)\ Sv (t) 



SP xy (t) = 2P xy ,-^- - - nmvlX . (16) 

v o,s * \ aa / a * =a * ^o,s 

It is convenient to introduce a generalized shear viscosity rj*(a*) through 

PUa*) = -V*(a*)a*. (17) 

In the Navier-Stokes (NS) approximation, 77* does not depend on a* by definition. Its 
dependence outside that limit is referred to as the rheological effects on the viscosity of the 
system. Use of Eq. (117]) into Eq. (TIB]) yields 



1 2 , 2 f dy*(a* 
D xy ,s+ -nmv QtS a s [—^r 



a*=a* 



where it has been used that 

6v (t) _ 1 ST(t) 
vo, s 2 T s 



(19) 



Next, the second term on the right hand side of Eq. ([131) will be considered. Substitution of 
the scaled form of the normal solution, Eq. ([7]) into Eq. (jSJ) yields 



with the dimensionless cooling rate (* given by 

(1 — a 2 )^^ 



(20) 



C (a 



2T (^) d 



dci I dc 2 c 12 x(c u a*)x(c 2 ,a* 



Therefore, 



sat) 



VQ,s 

2A 



c«)-a: 



d(*(a* 
da* 



8T(t) 



(21) 



(22) 



Now Eqs. (fl8l) and (1221) will be substitute into Eq. (fl3l) . A dimensionless representation will 
be employed. The relative deviation of the temperature 

5T(t) 



0(t) 



(23) 



is introduced, as well as a dimensionless time scale s defined by 

A 



(24) 



The time s is proportional to the average number of collisions per particle in the time interval 
between and t. In this way, Eq. ffT3"j) takes the form 



d6{s) 
ds 



C(a*) + £ 



a*=a* 



2 1 da* 



s . 



(25) 



Upon deriving this equation, the steady relationship given in Eq. ()3]) has been taken into 
account. It follows that a small homogeneous perturbation of the temperature of a granular 
gas under steady uniform shear flow decays exponentially in time, 



9(s) = e" 7 ^(0), 



(26) 



with the relaxation rate 



, *3 



7 = CK) + ^- 



dr]*(a* 
da* 



2 1 da* 



(27) 



a*=ai 



Note that, because of the constant linear relation between the dimensionless time scale s 



and the original one t, the decay of the perturbation in the latter scale is also exponential. 



III. MOLECULAR DYNAMICS SIMULATIONS 



Equation ( 126|) follows directly, without any additional approximation, from the hypothesis 
that the linear relaxation of the system after the perturbation is described by a normal 
distribution having the scaling form ([7]). But, to actually compute the decay rate 7, the 
dependence of both the cooling rate and the shear viscosity on the shear rate are needed. 
This would require to know the explicit form of the scaled distribution x(c, a*). In the NS 
approximation, the shear viscosit y is independent of the shear rate as already mentioned 



and, for a dilute granular gas 
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16j . the dependence of the cooling rate on the shear rate 



has been shown to be negligible in the linear approximation. Therefore, to NS order the 
rate 7 can be accurately approximated by the reduced cooling rate Q computed to zeroth 
order in the gradients, i.e. the one of the homogeneous cooling state particularized for the 
densit y a nd temperature of the steady USF state. An accurate expression for this property 



is 
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Co r (d/2) d v is ; ' [ } 

_ 16(1 - q)(l - 2a 2 ) 

C2{a} ~ 9 + 24d + (8d - 41)a + 30a 2 - 30a 3 ' ( } 

In the steady USF state, the inelasticity, measured through the coefficient of normal resti- 
tution a, and the shear rate a* are not independent, but they are related by Eq. (J3]). As 
a consequence, the restriction to small gradients required by the NS approximation, also 
implies limitation to small inelasticity, or to the quasielastic limit a — > 1. The analysis was 
consistently extended to the next order in a, Burnett order, in ref. joj. The general case 
of arbitrary shear rate (and inelasticity), has been studied by considering a simple kinetic 
model in which the inelastic Boltzmann collision operator is replaced by a single relaxation 
towards local equilibrium term, plus another term describing the energy dissipation due to 
the inelasticity of collisions (5), uj^ . Also, the Grad approximation to the Boltzmann equa- 
tion has been used [19j. Quite interestingly, all the above studies lead to expressions for 
the generalized shear viscosity that are hardly dis ting uishable over the complete interval 
< a < 1 and that are in very good agreement j^, Qjwith simulation results obtained by 
means of the direct Monte Carlo simulation method [2jJ . Moreover, another general feature 
of all the mentioned previous studies is that they do not lead to any shear dependence of 
the cooling rate of the steady USF state. Although we are not aware of any direct test of 
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this theoretical prediction by means of numerical simulations, the agreement observed for 
the steady temperature provides a strong support for it, since the expression of T s involves 
both the shear viscosity and the cooling rate (see Eq. ([3])). 

An accurate expression for the generalized shear viscosity of the time-dependent uniform 
shear flow state has been obtained in refs. Q]- In our notation, it reads, 

(d + 2)r(d/2) 



77V) 



2V2ir^/3 [1 + 2G(a*)} 2 [2 + H(a*)} 



where 



G(a*) = - sinh 2 < - cosh 1 
3 6 



1 + 



27 ( (d + 2)T{d/2)a* 



H(a*) 



(d+2)(l-a 2 ) 
4/3d 



32dn d - 1 V 
- 2G(a*)] [1 - 6G{a*)} 



[l + 6G(a*)Y 
1 + a 



(30) 



(31) 



(32) 



(33) 



For the sake of completeness, also the expression of the steady reducedd shear rate is written 
down, 



4^271^ 
(d + 2)T(d/2) 



d(d + 2)(l-a 2 ) 

8fid 



1/2 r 



a 



(34) 



Finally, d(*(a*)/da* will be also neglected in Eq. ( 1271) . Again, this is consistent with the 
previous model studies, but no numerical test of it has been carried out up to now, to the 
best of our knowledge. In summary, we approximate the expression for the rate of the decay 
of the temperature perturbation in Eq. ( l2Tj) by 

'a* 3 dr]*(a*) 



7 ~ Co + 



d da* 



(35) 



a*=a* 



with Co given by Eq. (T5H]) and the shear viscosity r/*(a*) by Eq. ([30]). For d = 3, Eq. (j35]) is 



the same as obtained by particularizing the expressions derived in ref. [22| by means of an 
extension of the Chapman- Enskog expansion to inhomogeneous steady states [23| , applied 
to the simple relaxation kinetic model equation mentioned above. 

In order to test the accuracy of Eqs. (1261) and (1351) . we have carried out event driven 
molecular dynamics (MD) simulations |24| of a system of iV = 2000 hard disks using Lees- 
Edwards periodic boundary conditions 25J. A square box was used and the density of the 
system was n = 0.02cr -2 , which is low enough as to expect the low density limit analysis to be 
accurate. The value of the shear rate was in all cases a = 6.32 x lO" 3 (T(0) /m) 1/2 tr -1 , where 



T(0) is the initial temperature. Finally, the value of the coefficient of normal restitution was 
varied in the simulations. 

Once the steady USF state was reached in the simulations, the small temperature pertur- 
bation was implemented as follows. The y-component of the velocity of all the particles was 
multiplied by the constant factor 1.2, i.e. the velocities were modified as v y — > 1.2v y . Note 
that this scaling of the velocity does not modify the macroscopic velocity field, keeping in 
particular a vanishing average value of its y-component. Therefore, this perturbation belongs 
to the class of perturbations considered in Sec. [Til In the simulations, an exponential decay 
of the temperature after the perturbation was observed, accordingly with the theoretical 
result in Eq. ( I26p . The single exponential shape of the decay guarantees that underlaying 
equation for the relaxation is linear and hence the validity of the linear approximation. In 
addition, it has been checked that equivalent results are obtained for smaller values of the 
perturbation of the velocity. Then, by means of a numerical fitting, the decay rate 7 was 
measured. The comparison of the simulation results and the theoretical prediction given in 
Eq. ( )35l) is shown in Fig. [TJ A fairly good agreement is observed except for very strong 
inelasticity (a very small). The discrepancy is probably due to the system approaching the 



instability point of the steady state 26] . To stress the relevance of the rheological effects on 
the rate, the Navier-Stokes rate, given by the dimensionless cooling rate Q, has also been 
plotted. 

An alternative way of measuring the decay rate 7 follows if the Onsager hypothesis on the 
regression of the equilibrium fluctuations |7| is extended and assumed to apply also to the 
internal energy of the steady USF state of a granular gas. Then, the decay of the internal 



energy time-correlation function in the steady state is also 7 |27j. The values obtained in 
this way have also been included in Fig. [TJ They agree with the results got from the decay 
of the initial perturbation, within the statistical uncertainties, with some discrepancy in the 
very inelastic limit. Again, this can be due to the proximity of the instability of the steady 
state, which can affect in a different way the time fluctuations and the linear response. 

IV. FINAL COMMENTS 

In this paper, the decay of an infinitesimal homogeneous perturbation of the temperature 
of a dilute granular gas in the steady USF state has been investigated, by using scaling 
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FIG. 1: (Color online) Dimensionless decay rate 7 of an small homogeneous perturbation of the 
temperature of a dilute granular in the steady uniform shear flow state as a function of the coefficient 
of normal restitution a. The solid line is the theoretical prediction given by Eq. (|35p . the dashed 
line is the dimensionless reduced shear rate Q, the (black) squares are the MD simulation results 
from the decay of an initial perturbation, and the (red) circles have been obtained from the decay 
of the energy time-correlation function in the steady state. 

properties of the normal solution of the Boltzmann equation. It has been shown that relevant 
rheological effects show up beyond the quasi-elastic limit. They tend to slow down the 
increase of the relaxation rate as the restitution coefficient decreases. Moreover, the rheology 
in the energy decay is associated to the shear rate dependence of the shear viscosity, while 
the dependence of the cooling rate, if it exits, it plays a negligible role. This confirms the 
results derived previously on the basis of model kinetic equations. 

Also, it has been verified that the energy time-correlation function in the steady USF 
state decays with the same rate as the energy perturbation, then obeying a hydrodynamic 
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description. This generalizes to a far from trivial situation the Osager hypothesis on the 
regression of equilibrium fluctuations. 

The natural question to be addressed in the near future is whether the conclusions reached 
here also apply to global properties of other systems in non-equilibrium steady states, both 
in granular gases and molecular fluids. If this were the case, it would indicate the existence 
of some generalization to non-equilibrium of the Onsager principle. This would be consistent 
with the structure of the fluctuating Navier-Stokes equations recently derived for dilute gran- 



ular gases [28] . Also, a few related studies concerning the validity of fluctuation-dissipation 



relations in driven dilute granular gases have already been carried out 
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Hail. 
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